A Rotating- Valence-Bond scenario for the 2D Antiferromagnetic Heisenberg Model 
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We propose that the valence bonds forming the ground state of the 2D-AF Heisenberg model on a 
square lattice may rotate under the effect of the antiferromagnetic background. To test this idea 
we apply a real space renormalization group approach to construct a rotating-valence-bond ground 
state anstaz. Our results are analytic and valid for any value of the spin S of the AF-magnet, which 
allow us to perfom semiclassical expansions of the energy and the staggered magnetization. 
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Despite of some initial controversies there is by now sumcient theoretical and experimental evidence, for the existencethe previou 
of antiferromagnetic long range order (AF LRO) in the 2d spin 1/2 Heisenberg antiferromagnetw (and -referenceshas total sp 
therein). This property has been observed in parent compounds of hight-T c materials such as La-JJuO^. From aspin 3S afte 
\ theoretical point of view this means that the strong quantum fluctuations implied by the low dimensionality andof bonds be 
■ the spin 1/2 do not destroy completely the Neel order, as it happensfl in Id. Though there is no a satisfactorysee below, r 
physical explanation of this fact, which may be important regarding the interplay bgtsveen antiferromagnetism andof long ranj 
t-H , superconductivity upon doping. The RVB scenario originally proposed by Andersonoa, while yielding an appealing To study 
picture of the ground state, does not explain the presence of AF LRO. This type of order may however be incorporatedFig.[l] a) as 
a posteriori in long range RVB ansatzs of factorized forrrjj, with gredictions similar to the ones obtained using Quantumlattice spac 
Monte Carlo methodsQ and variational plus Lanczos techniques!! A class of-,ghysical systems where the RVB approachby an effect 
may be actually realized is in spin ladders with an even number of chainsuu. The previous works leave still room to3 sites, obt; 
investigate in more depth the interplay between the RVB scenario, or more generally "valence bond scenarios" , andln our case, 
the AF order present in the 2d AF-magnets, described by the AF Hamiltonian H = JJ2{ij) &i ' to e = — J, 

In this letter we shall study a new possible effect of the antiferromagnetism on the dynamics of bonds, namely bond-perturbatio 
rotation. Let us recall that the basic mechanism considered by the RVB pictup. is the resonance between parallel The reno 

S neighbour bonds, which causes a substantial reduction of the ground state energylij. In addition to this effect we shallg^jg^ ^ j 
explore the possibility that the bonds also rotate around their ends under the influence of the AF background. To testgg m terms 
, this idea we propose a variational ground state in which the bonds rotate but do not resonate among themselves. This sum ru ^ e 
simplified ansatz allows us to decouple the rotation of bonds of other phenomena and in particular the resonance. A 
more realistic ansatz would include both rotation and resonance of bonds, but we shall not deal with this issue here. 
O ■ The analysis of the pure rotational-valence-bond ansatz leads us to a value of the staggered magnetization closer to 
the exact numerical result of the AFH model than to the Neel magnetization while we find the opposite result for 
the ground state energy. We interpret these results in terms of the previous picture of a bond scenario compatible 
with AF LRO and not containing the resonating mechanism, which is responsable for the lowering of the energy. An 
important ingredient of our construction is the use of real space renormalization group techniques, which allows us 
to obtain exact analytical results for any value of the spin S of the model (S is integer or half-integer and in the w here if a 
discussion above S = 1/2). The physical reason for this is that every cluster of 5 spins which share a rotating bond, the 5-block 
behave as an effective spin 3S 1 , coupled with its neighbours through a Heisenberg interaction with an effective couplingAmazingly 
constant. Hence the effective spin renormalizes to infinity, which allows us to compare the staggered magnetization 
and energy of the ansatz with the semiclassical spin- wave 1/S expansions of these quantities for the ground state of 
the AF-magnet. This comparison confirms the above picture of the rotating-valence-bond state as a state close to the 
exact ground state of the AF Heisenberg model in staggered magnetization but not in energy. 

Let us begin our approach by considering the cluster of 5 spins 1/2 of Fig.^ a). The configuration showed in Fig.^p or S pj n g 
a) is the exact ground state of the Ising piece of the Heisenberg Hamiltonian, given by H z = JY2i=i 4 ^o^i > w ^ ere (recall S — 
Sq and S 1 ? are the third component of the spin operators at the center and the i th position off the center respectively, ^he rq_ 
As soon as the "transverse" Hamiltonian H xy = JJ2i=i 4 $ x + S^Sf is switched on, the down-spin in the middle//' w hich t 
starts to move around the cluster, and a valence bond between the center and the remaining sites is formed in aj|Q_ e q Ua y 
s-wave (I = 0) symmetric state as shown in Fig.[I] b). Other rotational states with I ^ may appear corresponding to 
excitations (I being the orbital angula momentum of the bond). An alternative description of this state is given by 
first combining the 4 spins sourrounding the center into a spin 2 irrep, which in turn is combined with the spin 1/2 at 
the center yielding a spin 3/2 irrep with energy eo = — 3J/2. If instead of the spin 1/2 at each site there is a spin S 
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N' = ^, S' = 35, J' = 3£ 2 (5)J (3b) The abw 
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where the first contribution in Eq. ( |3a| ) comes from the energy of the blocks. As 3£ 2 (5) < 1, the flow equation ([3b]) o^AFH 6 irn 
implies that the coupling constant flows to zero J(") which means that the AFH model remains massless for a ^ ow f or re , 
arbitrary value of the spin 5. This fact allows us to compute the density of energy 600(5) (per site) as the following ne jgj 1 | :)0ur j. 
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Using eqs. (4) and (6) we can compute the ground state energy of our variational RG state for any value of the spin 

S. In particular for S=l/2 we get the value = -0.5464. This value has to he compared with the "exact" numerical 

result -0.6692, which is obtained using Green-function Monte Carlo methodstl We observe two facts: i) the error of 

the computation is 18%, which is surprisingly bigger than in Id whereto it amounts to a 12% and ii) the result is closer 

to the Neel state energy than to the exact ground state energy of the AFH model. This last observation is confirmed 

by the semiclassical expansion of eoo in the suhi S, which we compare with the standard formulas of Anderson and * Q n i eave , 

Kubo obtained using the spin wave methodsLj, 1 e. Manou 
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e s ^ = -25(5 + 0.158 + » + •••) 4 S.A. Kive 
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Observe that the lowest order correction to the Neel energy is absent in our case. 6 J- Carlsor 

Yet another explanation for the difference between our energy and the numerical value is that the rotating-valence- Hebb and 
bond state as depicted in Fig]2]does not have parallel adjacent bonds, which excludes the resonance among themselves. °- White, 
In order to have a better insight into the physics of the model it is convenient to compute the staggered magnetization^ ^- Dagott 
M = (jj ^2j(—^ySj). We have been able to obtain a closed formula for arbitrary spin 5 which is capable of analytical n Fazeka; 
study. To this purpose, we use the RG-equality for V.E.V. (^olCl^o) = {-0 O JC l^o) f° r renormalized observables O' ^riesrtcT) 
in the ground state and divide the sum in M into 5-block contributions. With the help of the renormalization spini 2 T n 
factors we arrive at the RG-equation for the staggered magnetization, L t r N 
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The explicit knowledge of £(5) (|l|) allows us to solve this RG-equation for the staggered magnetization in the ther- P arr i n 
modynamic limit N — > 00. In fact, as we know by now that the Hamiltonian renormalizes to its classical limit, we 
have lims^oo M(S) = S. Defining M(S) = 5/(5), Eq. (|) amounts to solving the equation /(5) = / (35) nf ™ 

subject to the boundary condition /(oo) = 1. Thus, we obtain the following formula for the staggered magnetization 
for arbitrary spin, 
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This is a nice formula in several regards. Forspin 5 — ^ we get ) — 0.373 to be compared with CL34i0.01 btained ro ^ a ^j n __ va jg 
with Green-function Monte Carlo methodso and Variational Monte Carlo plus Lanczos algorithm^. It amounts to a 

such as, e.g., spin wave theory plus 



7% error. Other approximate methods employed so far lead to values of M(^_ 
1/5-expaBsion givesLj 0.303, spin wave theory plus perturbation theory givesEE 
one foundt3 with pertubation theory around the Ising model to order 4 which is 0.371. Another interesting feature of 
our formula (Q) is that it allows us to make a 1/5-expansion yielding the result, 



0.313, etc. Our value is close to the 



M(S) = 5-0.2 + 0.06^ + <3(l/5 2 ) (8) 

which is to be compared with the spin wave theory result M (5) = 5 — 0.198 + 0(l/5 2 ) showing excellent agreement 
for the first two terms while discrepancies start from the term 1/5 onwards. 



2 



